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FOREWORD 


In  this  report,  the  authors  describe  an  efficient  computer  program  to  compute  geocentric 
tides  along  satellite  tracks  from  prepared  harmonic  tidal  constants  computed  on  a  standard  satellite 
track  grid  system.  The  program  is  an  improved  version  of  a  preliminary  program,  which  has  been 
applied  to  compute  instantaneous  geocentric  tides  along  SEASAT  tracks. 
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sponsorship  of  Mr.  O.  F.  Braxton,  Head  of  the  Strategic  Systems  Department,  his  Associate, 
Dr.  R.  J.  Anderle,  and  Mr.  C.  W.  Duke,  Jr.,  Head  of  the  Space  and  Surface  Systems  Division. 
Many  critical  and  stimulating  suggestions  were  gratefully  received  from  the  authors’  colleagues, 
Drs.  C.  J.  Cohen,  and  B.  Zondek. 
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ABSTRACT 


In  the  following  report,  the  authors  present  three  computer  programs  to  compute  geocentric 
tides  including  ocean  and  earth  tides  and  ocean-loading  effects  along  “parallel”  satellite  tracks 
from  the  harmonic  ocean  tidal  constants  listed  on  the  NSWC  GOTD  Tape  (Schwiderski  and 
Szeto  1981).  The  first  program  prepares  a  basic  standard  satellite  track  (SST),  which  is  shifted 
parallel  to  itself  to  define  a  1°  x  1°  SST  grid  system.  The  second  program  computes  harmonic 
geocentric  tide  data  at  all  SST  grid  points,  which  are  stored  on  three  magnetic  tapes  (SST  GTD 
Tapes  I,  II,  and  III).  From  the  prepared  SST  GTD  tapes,  the  third  (STT)  program  computes  in¬ 
stantaneous  tides  along  parallel  satellite  tracks  by  a  smooth  and  fast  interpolation  scheme.  All 
three  programs  eliminate  various  input-error  possibilities  of  their  preliminary  versions,  which 
have  been  applied  to  compute  instantaneous  geocentric  tides  along  SEASAT  tracks.  The  program 
descriptions  include  corresponding  User’s  Guides  and  Program  Listings.  An  extended  version  of 
the  STT  program  is  in  preparation.  It  will  include  group  beat  effects  on  all  major  tidal  components 
by  frequency-wise  neighboring  minor  tidal  modes. 
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I.  INTRODUCTION 


In  an  earlier  report,  the  authors  (Schwiderski  and  Szeto  1981)  described  briefly  the 
major  features  of  the  NSWC  ocean  tide  models  (Schwiderski  1978a,  b,  1979a  -  e,  1980, 
1981a  -  k,  1982a  -  d),  which  require  indispensable  consideration  in  various  applications.  In 
particular,  suggestions  were  discussed  to  improve  the  model  accuracy  especially  in  coastal 
waters.  These  general  discussions  were  followed  by  a  detailed  description  of  the  arrange¬ 
ment  and  format  of  the  NSWC  Global  Ocean  Tide  Data  (GOTD)  Tape  (Schwiderski  1981k), 
which  contains  the  1°  x  1°  gridded  harmonic  ocean  tide  constants,  i.e.,  amplitudes  and 
phases. 


Using  the  NSWC  GOTD  1981  tape,  a  Random-Point  Tide  (RPTIDE)  program 
was  elaborated  complete  with  User’s  Guide  and  Program  Listing.  The  RPTIDE  program 
computes  oceanic  and/or  geocentric  (including  earth  tides  and  loading  effects)  tides  at 
randomly  specified  geographical  points  and  instances.  Since  the  required  GOTD  1981  tape 
contains  over  one  million  data,  the  general  RPTIDE  program  is  cost-wise  limited  to  a 
relatively  small  number  of  random  points.  It  is  definitely  far  too  expensive  and  time 
consuming  to  compute,  for  instance,  geocentric  tidal  heights  along  satellite  tracks,  which 
carry  altimeters  to  measure  the  instantaneous  sea  surface  height  underneath  the  satellite 
every  half  second  or  so. 

It  is  the  purpose  of  the  present  report  to  present  the  special  NSWC  Satellite  Track 
Tide  (STT)  program,  which  computes  efficently  geocentric  tidal  heights  along  tracks  of 
altimeter-carrying  satellites  at  the  instances  sea-surface  measurements  are  being  taken.  It 
makes  effective  use  of  the  fact  that  the  ground  tracks  of  satellites  with  identical  and 
fixed-orbit  parameters  are  essentially  “parallel”  to  each  other  (see  Figure  2).  Indeed, 
disregarding  negligible  deviations,  two  consecutive  satellite  orbit  tracks  are  congruent  to 
each  other  in  space  and  time,  they  differ  only  in  a  uniform  longitudinal  westward  shift. 

In  Section  2,  one  finds  a  detailed  description  of  the  NSWC  SST  program  and  grid 
system,  which  generates  from  an  approximately  given  standard  satellite  orbital  ground  track, 
a  (basic)  SST  upon  which  the  basic  SST  grid  system  is  defined.  Section  3  presents,  in 
detail,  the  NSWC  Standard  Satellite  Track  Geocentric  Tide  Data  (SST  GTD)  program, 
which  converts  the  NSWC  Global  Ocean  Tide  Data  (GOTD  1981)  from  its  1°  x  1°  spheri¬ 
cal  grid  to  the  new  SST  grid  system  defined  in  Section  2.  The  generated  new  geocentric 
tide  data  are  stored  on  three  magnetic  tapes  (NSWC  SST  GTD  Tapes  I,  II,  and  III) 
described  in  Section  4.  These  tapes  are  used  in  the  NSWC  STT  program  described  in 
Section  5,  which  computes  instantaneous  geocentric  tides  along  specified  congruent  satellite 
tracks.  Finally,  the  Appendixes  A,  B,  and  C  contain  the  corresponding  User’s  Guides  and 
Program  Listings  of  all  three  NSWC  programs  (SST,  SST  GTD,  and  STT  programs). 
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2.  NSWC  STANDARD  SATELLITE  TRACK  (SST)  PROGRAM  AND  GRID  SYSTEM 


The  NSWC  SST  program  generates  from  a  given  approximate  standard  satellite 
ground  track  between  two  consecutive  ascending  nodal  points  a  basic  SST,  which  serves  as 
the  basis  of  the  SST  grid  system  defined  below. 


Definition  of  ESST  and  ASST 

An  “exact”  standard  satellite  track  (ESST)  is  defined  as  a  nonequatorial  one- 
revolution  ground  track  of  a  satellite  (say,  SEASAT  or  GEOS-3)  traveling  westward  around 
the  earth  between  two  consecutive  ascending  nodal  points  both  “exactly”  on  the  equator. 
If  one  or  both  nodal  points  deviate  slightly  off  the  equator,  the  track  is  called  “approxi¬ 
mate”  standard  satellite  track  (ASST). 

A.  Input  Data  (ASST) 

An  ASST  is  specified  by: 

(1)  (Xj,  0j)  =  longitudes  (East,  0°  <  X.  <  360°)  and  latitudes  (North,  -88° 

<  ^  <  +  88°)  of  j  =  I,  2,  3,...,  381  ASST  points,  which  are 
uniformly  spaced  by  the  constant  travel  time  A r  =  f/380,  where 

(2)  4  =  orbital  period  (in  sec) 

The  input  data  (1)  and  (2)  must  satisfy  the  following  ASST  accuracy  conditions. 

(a)  The  X^,  ^  values  must  be  given  to  0.5  •  10" 3  degrees  and  seconds,  re¬ 

spectively. 

(b)  The  exact  ESST  condition 

01  ~  0381  =  0  (1) 

is  in  practical  applications  usually  not  fulfilled.  However,  the  mbtimum  approximate  con¬ 
dition 


i  |  <  0.5°  and  |03#I|<O.5O  (2) 

can  and  must  be  enforced  to  avoid  significant  losses  in  accuracy.  The  SST  program  checks 
this  condition  (see  B.  below)  and  rejects  the  given  ASST  in  case  of  its  violation. 

Note  1:  The  Pole  regions  are  excluded  (14, 1  <  88°)  to  avoid,  singularities  in  the  grid 
system.  The  381  spacing  points  along  the  ASST  have  been  chosen  to  make  the  geo¬ 
graphical  distance  between  two  consecutive  points  about  equal  to  a  one-degree  equatorial 
distance.  Of  course,  any  other  number  of  points  could  be  chosen  in  principle. 
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B.  Main  Computation: 

The  following  procedure  generates  from  the  given  ASST  an  ESST,  which  is  needed 
for  the  largest-integer  arithmetic  in  the  STT  program  of  Section  5.  This  procedure  simply 
shifts  all  track  points  along  the  corresponding  parabolic  tangents,  in  order  to  enforce 
Equation  (1)  while  maintaining  a  constant  but  slightly  adjusted  travel  time  between  the 
new  consecutive  track  points  (see  Figure  1).  At  the  same  time  the  program  shifts  the 
resulting  ESST  to  the  Equator-Greenwich-Meridian  intersection  and  augments  the  track  by 
one  skew-symetric  additional  point  at  each  end.  This  completes  the  desired  generation  of 
the  (basic)  SST  upon  which  the  definition  of  the  SST  grid  system  is  based  (see  Figure  2). 
The  added  two  points  simplify  the  practical  application  of  the  STT  program  of  Section  5. 

(1)  Check  for 

I  0,  I  ^  0.5  and  I  03  g  ,  I  ^  0.5 
if  violated,  reject  given  ASST,  otherwise  compute: 

At  =  r/380 

AX,  =  X,  -  X2  (+360  if  <  0) 

A01  =  02  "  01 

AX3  81  =  ^3  8  0  ~  ^3  8  1  (+360  if  ^  0) 

A03  8  1  =  03  8  1  ""  03  8  0 
Atx  «  Ar0 i/A0j 
At38i  55  At03  g  1  /A03  8  1 
9  -  P  -  (A t3  8 1  -  At,  ) 

At  =  (Ar3  g,  -  Ar  1  )/3  8  0 

1  r  th 

v  =  —  [(AX,)2  +(A0,)2J 
X  =  360-  X,  -  0,  AX,/A0,) 

a  1  ~  4  ^3  s  1  +  03  •  1 AX3  8 1  / A03  g  ij  mod  360,  (0  <  X*  <  360) 

(2)  For  j  88  2, 3, ... ,  380  compute  consectively: 

A Tj  =  AT|_j  +At',  (At,  see  above) 

AXj  «  Xj_,  -  Xj+1  (+360  if  < 0) 

A*j  =  *j+l  "  *j-i 
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VArj/KAXj)*  +  (A^)2]* 

=  (X.  +  S.  AXj  +  X)  mod  360,(0  <  X'  <  360) 

“  "  SJ  A*i 


1 

EQUATOR 


Figure  1.  Construction  of  ESST  from  ASST 


•  Approximate  SST  points 


©Exact  points 
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(3)  Now  generate  the  (basic)  SST 

X,  =  360  -  X3 ,  <pi  —  -  03 

X2  =  360,  0 2=0 

Xj  =  X'  ,  0 .  =  0?  0  =  3,  4,...,  381) 

X382  =  X382,  0382  =  0 

X383  =  2X382  _  X381,  0383  =  -  0381 > 

^  =  720  -  X382  > 

and 

Ar  =  r'/380 

(4)  Print  all  SST  points  (Xj,  0j)  j  =  1,2, ... ,  383,  X,  and  Ar  on  the  SST  C 
Tape  I  as  described  in  Section  4. 

Note  2:  The  User’s  Guide  and  Program  Listing  are  given  in  Appendix  A. 

Definition  of  1°  x  1°  SST  Grid  System 

A  1°  x  1°  SST  grid  system  of  the  earth  is  defined  by  shifting  the  basic  SST  (see 
(3)  above)  one-degree-wise  westward  “parallel”  to  itself  as  shown  in  Figure  2.  Accordingly, 
the  grid  points  are  defined  by: 

Xjk  =  Xj  +  1  -  k  (+  360  if  <  0)  j  =  1,2, ... ,  383 

0jk  =  0j  k  =  1,2,...,  360 

Note  3:  As  can  be  seen  in  Figure  2,  the  SST  grid  system  is  not  a  unique  coordinate 
system.  Indeed,  through  every  given  point  in  the  SST  range,  one  finds  two  crossing  SST 
(ordinate)  lines  and  in  the  overlapping  region  around  the  equator  even  three.  Nevertheless, 
since  only  points  on  shifted  SSTs  will  be  considered  in  the  STT  program  of  Section  5,  a 
unique  geographical  orientation  is  possible  with  the  help  of  the  travel  time  Ar  between 
grid  points  on  the  SST. 
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"PARALLEL"  SST'S 


k  =  2  k  =  1 


EQUATOR  <p  =  0 


k  -  3  k 


Figure  2.  1°  x  1°  SST  Grid  Scheme 


(•)  “Paraflel”  Points 


i 


NSWC  TR  81-264 


3.  NSWC  STANDARD  SATELLITE  TRACK  GEOCENTRIC 
TIDE  DATA  (SST  GTD)  PROGRAM 


The  NSWC  (SST  GTD)  program  interpolates  at  all  SST  grid  points  (Section  2)  amplitudes 
and  phases  of  any  desired  partial  ocean  tide  by  second-order  ruled  surfaces  (linear  in  both  east 
and  north  directions)  using  the  NSWC  GOTO  tape  described  in  Schwiderski  and  Szeto  (1981). 
The  interpolated  values  are  modified  by  harmonic  addition  of  the  corresponding  amplitudes  and 
phases  of  the  earth  tide  and  of  the  earth  dip  in  response  to  the  ocean  tidal  load.  These  modi¬ 
fications  are  used  in  the  form  of  simple  Love  and  Accad-Pekeris  approximations  as  explained  in 
Schwiderski  and  Szeto  (1981).  Hence,  the  generated  new  harmonic  constants  on  the  1°  x  1° 
SST  grid  system  constitute  amplitudes  and  phases  of  the  total  geocentric  partial  tide. 

Note  1:  No  special  effort  is  made  to  improve  the  accuracy  of  the  oceanic  tides  in  coastal 
waters.  If  higher  accuracies  are  desired  in  such  areas,  instantaneous  tidal  computations  should 
apply  local  refinements  as  suggested  in  Schwiderski  ( 1 98  lj)  and  Schwiderski  and  Szeto  (1981). 


A.  Input  Data 

(1)  (Aj,0j)=  longitudes  (East,  0°  <  X  ^  360°)  and  latitudes  (North,  -  89°  <  0.  < 

+  89°)  ofj=  1,2,3  ....  3&3  SST  points  generated  by  the  SST  program  of 
Section  2  and  printed,  e.g.,  on  the  NSWC  SST  GTD  Tape  I  described  in 
Section  4  (Pole  regions  excluded  !). 

(2)  ($lm  ,  8lm  )  -  ocean  tide  amplitudes  (in  m)  and  Greenwich  phases  (in  deg)  from  GOTD 

1981  tape,  where  m  =  1,2, ...  ,  360  (longitude  numbers)  n  =1,2,  ... ,  168 
(latitude  number),  and 

(3)  i  =  specified  mode  number  1  <  i  <  1 1. 

Note  2:  fc1  =  9.999,  51  =  999.9  on  land 

>m  ,n  m  ,n 

8'm  n  =  360°  =  0°  (phase  jump) 

(4)  Earth  tide  parameters  in  GOTD  Mode  Order  (see  Schwiderski  and  Szeto  1981): 
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Mj 

i=  1, 

II 

JO 

E,  =0.148  308 

s2 

i=2, 

•'a  =  2, 

Ej  =0.069  059 

K, 

i  =  3, 

v3  =  1, 

Es  =0.086  638 

o, 

i  =  4, 

v*  =  1, 

E4  =0.061  551 

n2 

i=5. 

v$  =2, 

E,  =0.028  396 

P. 

i  =  6. 

II 

E6  =0.028  668 

k2 

i  =  7. 

vn  =  2, 

E,  =0.018  791 

Q. 

i=8. 

v s  =  1. 

E,  =0.011  785 

Mf 

i  =  9, 

v9  =  0, 

E,  =  0.025  546 

Mm 

i  =  10, 

vt0  =  0. 

E,0  =0.013  480 

Ssa 

i-  11. 

II 

© 

E, ,  =  0.01 1  901 

Main  Computation 

( 1 )  T ransfer  to  core  memory 

(Ar0.)forj=  1,2,3...  383 
Set  j  =  1  and  k  =  1 
(2)  Compute 

X.  =  (X3  +  1.5  -  k)  mod  360,  (0  <  X  <  360) 
m  =  Int  [Xj]  +  1, 

J/  =  m  -  X , 

n  =  Int  [90.5  -  +  1, 

0=0 =n-  (90.5 
If  n  >  169  (Antarctica!),  set 


and  go  to  (3)  below. 

If  n  >  169,  transfer  to  core  memory 
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where 

m  -  1  =  0  -*■  360  and  for  n  =  1 69 
«!.-•,  16.  =4. 1S, -9.999 
“in-,.  =*L,  =999.9 

Check  for  land  points  and  replace: 

(a)  if{J„  =  9.999,  replace  0  -*•  0 

(b)  if  f *m  „  =  9.999,  replace  6  -*■  1 

(c)  if  (a)  and  (b)  hold,  replace  \p  -+  1 

(d)  if^m_l  n_,  =  9.999,  replace  6  -*■  0 

(e)  if  _  j  n  =  9.999,  replace  6  -*■  1 
(0  if  (d)  and  (e)  hold,  replace  •*  0 

If  (c)  and  (f)  hold  (Land!),  set 

=s;)k  =o 

and  go  to  (3)  below,  otherwise  interpolate  £.  k  on  the  ruled  second-order  surface  in  \j/  and  6 

Test  for  360°  phase  jumps  and  replace 


(•) 


if  6 


i 

m  ,  n  - 1 


-  8 

m ,  n 


if  6* 


m-  1 ,  n- 1 


8‘_ 


if  8*  .  -  6l  _ 

m-1,  n  m ,  n 


>  180,  replace  8lm ,  ,  -  «‘m ,  n_ ,  -  360 

<  -  180,  replace  8‘m #  ,  -  8*n  ,  +  360 

>  180, replace ~  360 

<-  180,  replace  6^ _ ,  n _ ,  +  360’ 

>  180,  replace  Sj,,.,  §  360 

<-  180,  replace  «Li,  n  ^  Sm-i,  n  +  360 
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Use  adjusted  8’s  to  interpolate  k  by  formula  (*)  with  £  replaced  by  8.  Now  compute  and  replace 
(ocean  loading  effect,  see  Schwiderski  and  Szeto  1981) 

,k-  0.9333  *jk 


Compute  earth  tide  amplitude  function 


Ej  cos2  for  vi  =  2 
Ej  sin2  <j>.  for  v.  =  1 

^  Ej  (3  cos2  0.  -  2)  for  v{  =  0 


Finally,  compute  geocentric  harmonic  constants  by  harmonic  addition  (cos  and  sin  arguments 
in  degrees). 

*  k.*>' + «j>‘ + cos  iv + >-, -  k  ♦  Di }  * 

,  ijjksin5jk-E  sin[K(X  -k+l)] 

8‘  =  tan-  1 - - - 

^ ,k  cos  8j >k  +  Ej  cos  li'j  (Xj  -  k  +  1)] 

where  0  <  8  <  27r  in  radians. 

(3)  If  j  <  383,  replace  j  *»  j  +  1  and  repeat  (2)  above. 

If  j  <  383,  print  all  data  (£{  v ,  8|  v ) 

j  =  1,2,...,  383  and  k  on  the  magnetic  tape  SST  GTD  I,  II,  or  III  depending  on  the  mode 
i  as  described  in  Section  4. 

If  k  <  360,  replace  k  -*■  k  +  1 ,  j  — ►  1,  and  repeat  (2)  above. 

If  k  >  360,  stop  program. 


Note  3:  The  User’s  Guide  and  Program  Listing  of  this  program  are  presented  in  Appendix  B. 
Evidently,  by  applying  this  program  for  all  mode  numbers  i  =  1,  2,  3, ... ,  11,  one  generates  the 
SST  GTD  Tapes  I,  II,  and  III  described  in  Section  4. 
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4.  NSWC  STANDARD  SATELLITE  TRACK  GEOCENTRIC  TIDE  DATA  (SST  GTD) 

TAPES  I,  n,  AND  HI 


The  NSWC  SST  GTD  Tapes  I,  II,  and  III  contain  the  geocentric  (including  earth  tide  and 
loading  effects)  harmonic  tidal  constants;  i.e.,  amplitudes  and  phases 


.  i  1,  2, ,  1 1 

kk-SJI  '*'-2 . 383 

'  I  k  =  1,2 . 360 


generated  by  the  SST  GTD  program  described  in  Section  3.  The  data  are  defined  on  a  1°  x  1° 
SST  grid  system,  which  is  defined  in  Section  2  on  the  basis  of  an  SST  specified  by  the  longitudes 
(East,  in  deg)  and  latitudes  (North,  in  deg),  respectively, 

»,,♦,)  0  =  1, 2,  ...,383)  (*•) 

which  are  also  listed  on  SST  GTD  Tape  1. 

Computationally,  the  tidal  data  (•)  are  arranged  by  modes  i(=  1,2,...,  1 1)  on  the  three  magnetic 
tapes  in  the  order  shown  in  Table  1  below. 


Table  1 :  Mode  Arrangement  on  SST  GTD  Tapes 


SST  GTD  I 

i  =  2  :  S2 

* 

II 

i  =  4  :  O, 

SST  GTD  II 

z 

II 

i  =  6  :  P, 

N 

X 

r- 

II 

i  =  8  :  Q, 

SST  GTD  III 

1  i  =  9  :  Mf 

i  =  10  :  Mm 

i  =  1 1  :  Ssa 

- 

In  each  mode  i(=  1,  2, ...  ,1 1),  the  tidal  constants  (*)  are  arranged  by  SST-numbers  k  (=  1 , 
2, ... ,  360)  in  consecutive  pairs  of  blocks  with  each  block  containing  384  words  j  (=  1,2, ...  ,  384). 
The  first  383  words  (in  Format  FI 0.8)  in  the  first  block  are  amplitudes  k  (in  m)  and  in  the 
second  block  Greenwich  phases  6j  k  (in  rad).  The  last  word  0  =  384)  in  each  block  (in 
Formate  I  10)  gives  the  SST  number  k  (=  1,2, ... ,  360)  of  the  block  pair. 

As  shown  in  Table  1,  the  SST  GTD  Tape  I  contains  two  additional  (final)  blocks  of  384 
words  j(=  1,  2, ... ,  384)  in  Format  F  10.6.  The  first  383  words  in  the  first  block  represent  the 
SST  latitudes  <(>■  (in  deg)  and  in  the  second  block  the  SST  longitudes  X.  (in  deg).  The  last  word 
j  =  384  in  the  first  block  gives  the  SST  spacing  time  At  (in  sec),  and  in  the  second  block  the 
periodic  longitude  shift  X  of  the  SST  (in  deg.). 
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All  data  have  been  blockwise  generated  on  the  corresponding  magnetic  tapes  (Table  1) 
by  the  BUFFER-OUT  Statement  on  the  CDC  6700  computer.  These  tapes  have  the  following 
standard  properties:  7  track,  BDC  form,  even  parity,  556  bpi,  and  unlabeled. 

Note:  On  land  (see  Schwiderski  1978c)  all  tide  data  (*)  are  set  to  zero  (see  Section  3),  i.e. 

£J  k  =  8]  k  =  0  for  land. 


NSWC  TR  81-264 


5.  NSWC  SATELLITE  TRACK  TIDE  (STT)  PROGRAM 


The  NSWC  STT  program  uses  the  SST  GTD  Tapes  I,  II,  and  III  described  in  Sections  3  and 

4  to  compute  efficiently  instantaneous  geocentric  tides  at  equidistant  points  and  instances  along 
“constrained  satellite  (say,  SEASAT)  tracks  CSTY’  that  are  essentially  parallel-displaced  segments 
of  the  basic  SST  defined  in  Section  2  (see  Figures  1,  2,  and  3).  The  computed  geocentric  tides 
include  the  ocean  tides  superposed  with  the  corresponding  earth  tides  and  earth  dips  in  response  to 
the  oceanic  tidal  loads  (see  Schwiderski  and  Szeto  1981).  When  all  leading  tidal  modes  (M2,  S2, 
K, ,  O, ,  N2,  P, ,  K2,  Q| ,  Mf,  Mm,  and  Ssa;  see  Table  2)  are  included,  the  resulting  instantaneous 
tidal  elevations  carry  a  10-cm  accuracy  anywhere  over  open  ocean  areas  (see  Schwiderski  1978a,  b, 
1979a-e,  1980,  1 98 1  a— j,  and  1982a-d).  This  accuracy  diminishes  somewhat  in  coastal  waters 
where  special  improved  computations  are  suggested  in  Schwiderski  ( 1 98 1  j)  and  Schwiderski  and 
Szeto  (1981). 

Definition  of  Constrained  Satellite  Tracks  (CST) 

A  CST  is  a  uniformly  spaced  (at  least  two  points)  segment  of  a  satellite  ground  track, 
which  (see  Figure  3): 

(a)  Is  almost  parallel  (congruent)  to  the  basic  SST  of  the  SST  GTD  Tapes  I,  II,  and  III 
defined  in  Sections  2,  3,  and  4. 

(b)  1  s  grid  wise  continuous,  i.e.,  it  is  without  gaps  in  equidistant  spacing  points. 

(c)  Lies  almost  entirely  between  the  track’s  two  consecutive  ascending  nodal  points 
on  the  equator. 

(d)  Lies  almost  entirely  over  the  global  ocean  area  specified  by  non-zero  tidal  constants 
listed  on  the  SST  GTD  Tapes  I,  II,  and  Ill  of  Section  4. 

Note  1:  If  condition  (a)  is  violated,  say,  by  more  than  0.5°  along  the  given  CST  when  shifted  to 
coincide  with  the  SST,  then  a  new  appropriate  SST  and  corresponding,  SST  GTD  tapes  must 
be  prepared  (see  Sections  2,  3,  and  4). 

Note  2:  Satellite  tracks  that  are  not  gridwise  continuous,  cross  the  equator,  say,  by  more  than 

5  sec  and/or  pass  over  land  areas  must  be  broken  up  into  separate  segments  of  CSTs  to  fulfill  the 
conditions  (b,  c,  and  d)  above. 

A.  Input  Data: 

(1)  y  =  year  >  1975  (fixed  for  one  run!) 

(2)  d  ■  day  of  year  y  (d  =  1  for  January  1st,  also  fixed  for  one  run!) 
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(3)  t(  =  initial  time  (in  sec)  of  first  CST  point  (Xt ,  )  relative  to  Greenwich  midnight 

of  day  d  (universal  time) 

(4)  At  =  constant  time  step  (in  sec)  along  CST 

(5)  N  ( >  1 )  =  total  number  of  CST  points 

(6)  (X„,  0„)  =  longitudes  (East)  and  latitudes  of  first  two  (n  =  1,  2)  CST  points  (in 
deg.) 

(7)  (Xg,  Xb )  =  equator-crossing  longitudes  (in  deg.)  corresponding  to  the  two  con¬ 
secutive  ascending  nodes  of  the  track  containing  the  CST 

(8)  (ta,  tb)  =  equator-crossing  times  (in  sec)  relative  to  Greenwich  midnight  of  day 
d  (d  +  1),  or  (d  -  1)  for  ta  belonging  to  Xg  and  of  day  d,  (d  -  1),  or  (d  +  1)  for  tb 
belonging  to  Xb 

(9)  (£!  k ,  5(  k )  =  geocentric  tidal  amplitudes  (in  m)  and  Greenwich  phases  (in  rad)  from 
S!>T  GTD  Tapes  1,  II,  and  III,  where 

i  =  1 ,  2, ...  ,  I  (<  1 1 )  =  mode  numbers 

j  =  1,2,  ...  ,  383  =  SST  spacing  points 

k  =  1,2, ...  ,  360  =  SST  ordinate  lines  (see  Figures  2  and  3) 

(10)  I  =  total  number  of  tidal  modes  ( 1  <  I  <  1 1 )  to  be  superposed 

(11)  A t  =  SST  time  step  from  SST  GTD  Tape  l  (in  sec) 

(12)  X  =  period  longitude  shift  (in  deg)  of  SST  from  SST  GTD  Tape  I 

Note  3:  The  CST  data  (3,  4,  5,  and  6)  must  satisfy  the  CST  definition  above,  necessary  splits 
will  be  requested! 

Note  4:  For  the  equator-crossing  longitudes  (Xa,  Xb)  and  the  corresponding  times  ta  and  tfc  error 
bounds  of  less  than  0. 1  °  and  1 .0  sec,  respectively,  are  strongly  recommended.  Errors  of  more 
than  0.4°  and  4.0  sec  will  be  rejected  for  corrections  (see  B(3)  and  (4)  below). 

Note  5:  The  tidal  data 

£■  k  =  k  =  0  signal  land,  and 

}‘  =  2n  =  0  a  phase  jump. 

j  i  * 
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Table  2.  Constants  of  Major  Tidal  Modes 


Tidal  Mode 

K  (m) 

a  (10  4  /sec) 

X (deg) 

Semidiurnal  Species 

M2  =  Principal  Lunar 

0.242  334 

1.405  19 

2h0  -  2sq 

S2  =  Principal  Solar 

0.112841 

1.454  44 

0 

Na  =  Elliptical  Lunar 

0.046  398 

1.378  80 

2h0  —  3s0  +  Po 

K2  =  Declination  Luni-Solar 

0.030  704 

1.458  42 

2hg 

Diurnal  Species 

K,  =  Declination  Luni-Solar 

0.141  565 

0.729  21 

h0  +90 

Ox  =  Principal  Lunar 

0.100  574 

0.675  98 

h0  -  2s0  -  90 

Pi  *  Principal  Solar 

0.046  843 

0.725  23 

-h0  -  90 

Qx  =  Elliptical  Lunar 

0.019  256 

0.649  59 

h0  -  3s0  +  p0  -  90 

Long-Period  Species 

Mf  =  Fortnightly  Lunar 

0.041  742 

0.053  234 

2s0 

Mm  =  Monthly  Lunar 

0.022  026 

0.026  392 

s0  -  Po 

Ssa  =  Semiannual  Solar 

0.019  446 

0.003  9821 

2h0 

K  -  amplitude  of  the  partial  tide 
o  =  frequency  of  the  partial  tide 
X  =  astronomical  argument  of  the  partial  tide 

(h0,  sQ,  p0)  =  mean  longitudes  of  sun,  moon,  and  lunai  perigee  at  Greenwich  midnight 
h0  =  279.696  68  +  36  000.  768  925  485T  +  3.03  •  10"4T2 
s0  =  270.434  358  +  481  267.883  141  37T- 0.001  133T2  +  1.9  •  I(T6  T3 
p0  =  334.329  653 +  4  069.034  032  957  5T- 0.010  325T2  -  1.2-  10“ 5  T3 

where 

T  =  [27  392.500  528  +  1 .000  000  035  6DJ/36  525 
D  =  d  +  365(y  -  1975)  +  Int  t(y  -  1973)/4] 
d  =  day  number  of  year  (d  =  1  for  January  1 ) 
y  >  1 975  =  year  number, 

and 


Int  [xj  =  integral  part  of  x 
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B.  Computation  of  SST  Grid  Brackets  and  Other  Constants  for  CST 

In  the  following  preliminary  computations,  the  given  CST  will  be  bracketed  between  two 
consecutive  SST  grid  ordinates  k  and  (k  +  1)  as  shown  in  Figure  3.  Subsequently,  time  and  space 
constants  will  be  computed,  which  are  needed  in  the  following  main  computations. 

(1)  Transfer  to  core  memory  the  CST  data 

y,  d,  tj,  At,  N,  ta,  tfe,  Xa,  Xb,  (Xj ,  ),  (X2, 02),  and  from  the  SST  GTD  Tape  1  the  SST  data 

A 

At  and  X. 

(2)  Compute  and  adjust  day  count 

X1  =  360  -  Xb  +  Xa  (+  360  if  <  360), 
t  =  ^  -  ^  (+  86400  if  <  0), 
t  =  380  At,  t.  =  t.  -  t  . 

1  J  a 
A 

lfr}  >  3t,  replace 

Tj  Tj  -  86400,  tj  tj  -  86400,  d  -+  d  +  1 . 

A 

IfTj  <-3t,  replace 

Tj  “►  Tj  +  86400,  tj  tj  +  86400,  d  -►  d  -  I. 

(3)  Compute  and  check 

rN  +(N-  DAt. 

If  Tj  <  -  At/2  and/or  tn  -  r  >  At/2,  stop  and  print:  Check  Track  Data!  Otherwise 
compute  and  replace 

At-*AtT/t,  t€  =  ta, 

A  A 

Tj  -►Tj  T/t 

j  A  (  +  360  if  <  0 
Xg  *  \  +  -  (X1  -  X)  <  -  360  if  >  360 
(  +  0  otherwise. 
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(4)  If  |t,  I  <  At  compute 

0  =  0j/(02  ~  )> 

AX,  =  X,  -  X2  (+  360  if  <  0), 

AX  =  Xe-  X,  (+360  if  <  -  10), 

AX'  =  0AX, , 
t'  =  </>  At. 

If  |t,  -  fj  |  >  0.4  At 

and/or 

| AX  -  AX ' '  >  0.4, 

stop  and  print:  Check  Track  Data!  Otherwise  replace 

te-te  +  T,  -t,', 

S+360  if  <  0 
-  360  if  >360 
+  0  otherwise. 

(5)  Compute  the  CST  bracket  data  (see  Figure  3) 
k  =  Int  {361  -  Xe] ,  X  =  (361  -  Xe)  -  k. 

Replace 

k  =  361  ^  1. 

Compute  and  replace 

J  =  lnt  [2  +  t, /At], 

M  =  Int  [4-  J  +  tn/AtJ, 

t,  -*■  t,  -  (J  -  2)At 
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(6)  Compute  the  constants 


At  =  1/At,  fj  =  At/2,  t2  -t,  At,t3  =  ?2  At 

With  the  tidal  parameters  hQ,  sQ,  p0,  and  (Oj,  Xp  i  =  1,2,...,  1 1)  listed  in  Tables  I  and  2  computed 
fori- 1,2,  ...,11. 

of  =  At  ,Xj  =  o.te-2  ot -MTXj/180 

Keep  the  constants  I,  J,  k,  M,  N,  X,  At,  Ar,  Ar,  t}  ,  r}  ,r^,r3,  and  (o. , x4 *  i  =  1,  2, ...  ,1 1)  for  the 
main  computations. 

C.  First  Rough  Interpolation 


In  this  first  step  tidal  amplitudes  and  phases  are  mode-wise  linearly  interpolated  at  the 
SST-spacing  points  on  the  CST  using  the  data  along  the  neighboring  SST  ordinate  lines  k  and 
(k  +  1)  as  shown  in  Figure  3.  Subsequently,  instantaneous  geocentric  tides  £  are  computed  from 
the  interpolated  harmonic  constants,  which  are  automatically  mode-wise  superposed, 

(1)  Set:  i  -  1  and  fm  «  0  for  m  =  1,  2,  ... ,  M 

(2)  Transfer  to  core  memory  the  tidal  constants 


(£j,k’  ^j,k  ^  an(*  1 9  ®j,k+I 


) 


forj  «  1,2, ... ,  383  and  i  and  k  fixed  (k  +  1  =  361  -►  1!) 


Set 


j  =  J,  m  =  1 
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(3) 


Check  for  land  points  and  define  ^  ,  X2 )  as  follows 
(  >  0  set  X2  =  X 


For  | 
FOr*j.k+i| 


>  0  set  X2  =  1 

>  0  set  X,  =  1  -  X 

>  0  set  X,  =  1 


Now  interpolate  linearly 


«.  =  M,k+M.k  +  . 

check  for  27r-phase  jumps  and  replace  6’s  for 

>7r  replace  6j  k+1  -*6jk+1  +  2ir 
<-  n  replace  6|  .  *-►  5j  .  +  2tt 
Use  adjusted  6’s  to  compute  6m  by  formula  (*)  with  £  -*  5 
Now  compute  the  superposed  instantaneous  tide 

jF  +  C0S  (ai  j+X;  -  Sm  ) 

If  m  <  M,  replace  m  -►  m  +  1 ,  j  -►  j  +  1,  and  repeat  (3)  above.  If  (m  <  M)>  go  to  (4)  below. 


5j,k  ~  ^j,k  +  1 


<*) 


(4)  If  i  <  I,  replace  i  i  +  1  and  repeat  (2)  above. 

If  (i  <  I),  follow  D  below,  but  keep  the  constants  M,  N,  At,  Ar,  A?,  r, ,  ,  t"2,  r^,  and 

fm  (m  =  1  >2, ...  ,  m). 

D.  Second  Refined  Interpolation 

The  finally  desired  instantaneous  geocentric  tides  £  (n  =  1,2,...  ,  N)  at  the  given  CST 
spacing  points  are  generally  computed  by  a  “cubic-parabolic  spline”  interpolation  to  achieve 
a  smooth  tangent  variation  along  the  fn.  In  order  to  compute  the  tide  f  ,  say,  between  the  SST  - 
spaced  data  and  fa  (see  C  and  Figure  3)  a  “cubic”  polynomial  is  forced  through  f2  and  f3 
with  the  corresponding  “parabolic”  slopes 

*m  =  + 1  ”  - 1  for  m  =  2,  and  3 
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Naturally,  for  shorter  CSTs;  e.g.,  between  land  areas  (signaled  by  ?  =  0  data)  only  ordi¬ 
nary  parabolic  or  linear  interpolation  is  used.  In  detail,  the  following  tests  must  be  made,  which 
lead  to  the  different  interpolation  cases  1,2,3,  4,  and  (subsequently)  5. 

( 1 )  If  M  >2  and  ?2  ^  0.  and 

(1.1)  ?3  ^  0,  and 

(1.1.1)  f^O,  set 

r  =  Tj ,  m  =  1 ,  n  =  1,  and  go  to  case  3 

(1.1.2)  (ft  =  0),  set 

r  =  r !  -  At,  m  =  2,  n  =  1 ,  and 

( 1 . 1 .2. 1 )  for  M  >  3  and  ?4  =£  0,  go  to  case  3 

(1.1. 2. 2)  otherwise  go  to  case  2 

(1.2)  (?3  =0  hand 

(1.2.1)  fj  *0,set 

r  =  m  =  1 ,  n  =  1  and  go  to  case  2 

( 1 .2.2)  (?!  =  0),  go  to  case  ( 1 ) 

(2)  If  (1)  fails  and 

(2.1 )  M  >  2,  ?2  ¥=  0,  am/?,  #  0,  set  r  =  r, ,  m  =  1 ,  n  *  1,  and  go  to  case  2 

(2.2)  otherwise  go  to  case  1 

Now  compute  the  interpolutions 

Case  1 .  Constant  interpolation  for  one  oceanic  datum  ?!  =£  0  or  ?2  =£  0 
Compute 

?n  =  ?!  +  f2  f°r  n  =  1 , 2, ... ,  N,  and  go  to  E 


i- 
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Case  2:  Linear  interpolation  for  two  oceanic  data  f m  #  0  and  fm  + ,  #  0 
With  a  =  &t~ (fm  +  1  -  fm  ),  b  =  a  At,  compute 


?  1  =  ? m  +  a  T-  tor  n  =  1 


fn  =  fn_j  +  b,  for  n  =  2,  3 . N,  and  go  to  E 


Case  3:  Parabolic  interpolation  for  three  oceanic  data  fm  ^  0.  fm  =£0,  andfm  +  2  0 

(a)  Compute: 

a  =  r7  ('3L  +  4fn  +  i  -  +  2  >'  b  =  ^  <fm  -  2?mtl  +  +2  >’  and  go  for 


M  <  m  +  3,  or  f  +3  =  0  to  (b),  otherwise  go  to  (c) 


(b)  Compute 


fn  =(br  +  a)T  +  fm 


If  n  <  N,  replace  n  n  +  1 ,  r  -*  r  +  At,  and  repeat  (b). 


If  (n  <  N),  go  to  E. 


(c)  Compute 


?n  ~  (br  +  a)  t  +  f 


Replace  n  -►  n  +  1  and  r  -►  r  +  At 


If  r  <Ar,  repeat  (c) 


If  (t  <  At),  replace  r  -►  r  -  At,  m  m  +  1 ,  and  go  to  Case  4  below . 

Case  4:  Cubic-parabolic  interpolation  for  four  oceanic  data,  =£  0,  =£  0,  fm  +  J  ^  0, 

andfm  +  2*0: 

(a)  Compute 

-?i  (fm.l  -  fm-,)<  b  =  ^  (2  fm-l  -  S  fm  +4fm  +  1  -  ?m+2) 

c  =  ?3  (fm  +  2  -  3fm  +  i  +3fm  -  rm„,).andgoto(b)below 


22 


NSWC  TR  81-264 


(b)  Compute 

fn  =  [(cr+b)r  +  a)j  r  +  fm 
Replace  n -*■  n  +  1  and t-*t  +  A t 
If  t  <  At,  repeat  (b) 

If  (t  <  At),  replace  r  -*■  r  -  At,  and  m  -*•  m  +  1 ,  and  if  M  <  m  +  2,  or  +  2  =0,  go  to 
Case  5,  otherwise  repeat  (a)  above. 

Case  5:  Parabolic  end-point  interpolation  for  three  oceanic  datafm_J  #  0,  fm  #  0,  and  ?m  +  ] 
#0: 

(a)  Compute 

a  =  *I  C  -  =  +  l  -  2fm  +  fm  _ , and  go  to  (b)  below 

(b)  Compute 

fn  =(br  +  a)r  +  fm 

If  n  <  N,  replace  n  -»•  n  +  1 ,  r  '-*•  r  +  At  and  repeat  (b) 

If  (n  <  N),  go  toE 
E.  Output  Data 

List  and/or  print  on  tape  all  tidal  data 

fn  :  n  =  1, 2, ,  N 

as  specifically  requested. 

Note  5:  The  User’s  Guide  and  Program  Listing  of  this  program  are  given  in  Appendix  C. 
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APPENDIX  A 

1.  SST  PROGRAM  USER’S  GUIDE 


The  Standard  Satellite  Track  (SST)  program  generates  from  a  given  Approximate 
Standard  Satellite  Track  (ASST)  (see  Section  2,  Part  A)  an  Exact  Standard  Satellite  Track 
(ESST)  that  is  used  by  the  SSTGTD  (see  Section  3)  program. 

The  following  listing  shows  the  SST  program  written  in  CDC  Fortran  Extended  for 
the  CDC  6700  computer  under  the  SCOPE  3.4  operating  system. 

Input  to  the  SST  program  consist  of  the  ASST  File  attached  to  TAPE1  and  a  data 
card  containing  the  orbital  period  in  seconds  with  E20.9  format.  The  computed  ESST  is 
written  onto  a  file  called  TAPE2.  This  ESST  will  be  called  the  SST. 
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APPENDIX  B 

1.  SSTGTD  PROGRAM  USER  S  GUIDE 


The  standard  Satellite  Track  Geocentric  Tide  Data  (SSTGTD)  program  generates 
geocentric  harmonic  constants  at  SST  (see  Section  2)  grid  points.  This  data  will  be  used 
by  the  STT  (Section  5)  program. 

The  following  program  listing  shows  SSTGTD  written  in  CDC  Fortran  Extended  for 
the  CDC  6700  computer  under  the  SCOPE  3.4  operating  system. 

Input  for  the  SSTGTD  program  consist  of  three  Files  ard  two  other  values  using  the 
free-format  READ  statement.  These  two  values  are  the  species  number  and  the  earth- 
tide  amplitude  of  the  tide  that  is  being  processed.  The  computed  values  are  put  on  two 
other  Files  (TAPE4  and  TAPE5). 

The  following  gives  the  control  cards  and  data  card  for  a  sample  run  of  the  SSTGTD 
program  for  the  M2  tide. 

(job  card  > 

(account  card) 

FTN,  R=3,  A. 

ATTACH,  TAPE2,  the  SST  file  generated  by  the  SST  (Section  2)  program. 

ATTACH,  TAPE1,  £‘m  n  (amplitude  in  meters)!  m  =  1,  2,...,  360  (longitude  number) 

>n  =  1,2,...,  168  (colatitude  number) 

ATTACH,  TAPE3,  8lm  (phase  in  degrees)  ]i  =  tidal  mode  being  processed 

Tapes  1  and  3  contain  GOTD  values 
made  randomly  accessible  by  colaiitude 
number  n. 

REQUEST.  TAPE4,  *PF. 

REQUEST,  TAPES,  *PF. 

LGO. 

CATALOG,  TAPE4,  R  (amplitudes  in  meters  K  j  =  1,2,,..,  383  (SST  points) 

CATALOG,  TAPE  5,  6j  k  (phases  in  radians)  jk  -  1,2,...,  360  (SST  number) 

i  =  tidal  mode  being  processed 

( end  of  record  > 

(SSTGTD  program) 

(end  of  record) 

(data  card)  for  i  =  1  (M2  tide); 


(2,  0.148308) 
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115  6050  FORMAT ( *  *,*SSTGTD  PROGRAM  COMPLETED*) 

CALL  EXIT 

12  STOP  -12. STOPPED  IN  SSTGTD“ 

13  STOP  -13. STOPPED  IN  SSTGTD" 
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N  INTEGER  F.P.  REFS  29  DEFINED  t 

PHI  REAL  F.P.  REFS  26  DEFINED  1 

P5I  REAL  F.P.  DEFINED  1  28 

P9  REAL  REFS  27  29  DEFINED 


NSWC  TR  81-264 


APPENDIX  C 

STT  PROGRAM  USER’S  GUIDE 
AND  PROGRAM  USTING 


i 

j 


NSWC  TR  81-264 


APPEN1BX  C 

1.  STT  PROGRAM  USER’S  GUIDE 


The  Satellite  Track  Tide  (STT)  program  computes  instantaneous  geocentric  tides  at 
equidistant  points  along  a  given  Constrained  Satellite  Track  (CST). 

The  following  program  listing  shows  the  STT  (program)  subroutine  written  in  CDC 
Fortran  Extended  for  the  CDC  6700  computer  under  the  SCOPE  3.4  operating  system. 

The  STT  subroutine  is  called  with  the  following  statement: 

CALL  STT  (MODES,  YEAR,  DAY,  TIME1,  NTOTAL,  IUNIT,  1SST, 

CSTLAT,  CSTLON,  ELONA,  ELONB,  ETIMEA,  ETIMEB, 

TIDE,  DT,  NPTS,  ISTOP) 


where 

MODES  =  the  number  of  tidal  modes  to  process 

YEAR  =  year  >  1975 

DAY  =  day  of  YEAR,  example:  February  1,  1978,  DAY  =  32 

TIME1  =  time  (in  sec)  of  1st  point  of  the  CST  to  be  processed  relative  to  Greenwich 

midnight  of  DAY 

NTOTAL  =  number  of  points  on  the  CST  to  be  processed 

IUNIT  =  beginning  unit  number  of  (2  x  MODES)  consecutive  units  to  which  the 
SSTGTD  files  are  attached 

ISST  =  unit  number  to  which  the  SST  File  is  attached 

CSTLAT  =  array  containing  latitudes  of  the  1st  2  CST  points  in  degrees 


CSTLON  =  array  containing  longitudes  (EAST)  of  the  1st  2  CST  points  in  degrees 
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ELONA  =|equator  crossing  longitudes  (in  deg)  corresponding  to  the  2 
ELONB  =  (consecutive  ascending  modes  of  the  track  containing  the  CST 

ETIMEA  =  )equator  crossing  times  (in  sec)  relative  to  Greenwich  midnight 
ETIMEB  =  lot  DAY,  (DAY-1),  or  (DAY+1)  for  HTIMEA  belonging  to  ELONA  and 
/ETIMEB  belonging  to  HLONB 

TIDE  =  array  containing  the  generated  tide  values 

DT  =  time  in  sec  between  consecutive  points  of  the  CST 

NPTS  =  number  of  tide  values  returned  by  STT  in  T1DE(1)  through  TIDE  (NPTS), 

1  <  NPTS  <  100 

ISTOP  =  flag,  if  ISTOP  =  0,  the  last  set  of  tide  values  have  been  generated  for  the 
given  CST. 

For  each  CST  track  to  be  processed,  the  STT  subroutine  (program)  is  called  repeated 
until  a  zero  is  returned  for  the  ISTOP  variable.  For  each  successive  return  from  STT, 
NPTS  indicates  the  number  of  successive  tide  values  generated  and  stored  in  the  TIDE 
array.  (1  <  NPTS  <  100)  Also,  NPTS  must  be  set  to  -1  before  the  first  call  to  STT 
for  each  track.  The  BUFFER  IN  statement  in  STT  bring  in  data  files  (SSTGTD)  that 
are  in  binary  form.  The  SSTGTD  magnetic  tapes  are  in  coded  form  as  described  in 
Section  4.  The  SST  data  file  is  also  described  in  Section  4.  The  following  shows  how  the 
SSTGTD  files  are  attached. 

ATTACH,  TAPE  8, 

ATTACH,  TAPE  8+1, 

ATTACH,  TAPE  8+2, 

ATTACH,  TAPE  8+3, 

* 

ATTACH,  TAPE  8+21, 

ATTACH,  TAPE  8+22, 


5i= 


=  1,2 .  11  (tidal  model) 

j  =  1,2 .  383  (SST  points) 

k  =  1,2,  360  (SST  number) 

SSTGTD  files,  see  Section  4 
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